Local load sharing fiber bundles with a lower cutoff of strength disorder 
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We study the failure properties of fiber bundles with a finite lower cutoff of the strength disorder 
varying the range of interaction between the limiting cases of completely global and completely local 
load sharing. Computer simulations revealed that at any range of load redistribution there exists 
a critical cutoff strength where the macroscopic response of the bundle becomes perfectly brittle, 
i.e. linearly elastic behavior is obtained up to global failure, which occurs catastrophically after 
the breaking of a small number of fibers. As an extension of recent mean field studies [Phys. Rev. 
Lett. 95, 12550f (2005)], we demonstrate that approaching the critical cutoff, the size distribution 
of bursts of breaking fibers shows a crossover to a universal power law form with an exponent 3/2 
independent of the range of interaction. 

PACS numbers: 46.50.+a, 62.20.Mk, 81.40.Np 
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In engineering constructions solids are subjected to 
various types of external loads and typically fail when 
the load exceeds a critical value. Monitoring stressed 
systems and forecasting an imminent failure event is of 
enormous importance due to the possible material and 
human costs. In spite of the large amount of experimental 
and theoretical efforts that has been undertaken over the 
past decades, there is no comprehensive understanding of 
failure phenomena, which is also reflected by the absence 
of reliable prediction methods. Materials of low disorder 
typically fail in a "one-crack" way, where the main prob- 
lem is to prevent crack initiation and propagation. How- 
ever, the failure of highly disordered materials proceeds 
in bursts of local breaking events which can be recorded 
in form of acoustic signals. Experiments on a large vari- 
ety of materials have revealed that in crackle noise spec- 
tra accompanying quasi-static fracture of disordered ma- 
terials, the amplitude and duration of signals and the 
waiting time between them are characterized by power 
law distributions over a broad range 0,1113. Quantita- 
tive changes of the burst activity when approaching the 
critical load could be precursors of catastrophic failure 
and may serve as the basis for forecasting techniques. 

In the framework of discrete models of the fracture of 
disordered materials, bursts can be identified as trails 
of correlated breakings of the microscopic constituents 
of the model. Fiber bundle models consist of a parallel 
bundle of fibers with identical linearly elastic behavior 
and randomly distributed breaking thresholds 0, IE El • 
Under an increasing external load, each fiber breaking 
is followed by a load redistribution over the remaining 
intact fibers, which may trigger avalanches of correlated 
breaking events analogous to crackling noise in experi- 
ments. Assuming global load sharing (GLS), it has been 
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shown that the distribution of avalanche sizes has a uni- 
versal power law behavior with an exponent 5/2 0- In 
the other extreme of local load sharing (LLS), redistribut- 
ing the load solely over the closest neighborhood of fibers, 
the avalanche distribution appears also to be a power law 
but with a higher exponent w9/2 0, • For global load 
sharing it has recently been pointed out that the distribu- 
tion of burst sizes significantly changes if the weak fibers 
are removed from the bundle: if the strength distribu- 
tion of fibers has a finite lower cutoff, or analogously, if 
the recording of avalanches starts after the breaking of 
the weak elements, the burst size distribution is found to 
show a crossover to another power law with a significantly 
lower exponent 3/2 0E3- The effective range of inter- 
action in real materials may have large variations 
therefore, in order to use the crossover effect of burst sizes 
in forecasting of imminent failure, its robustness with re- 
spect to the range of interaction has to be explored. 

In the present paper we extend recent mean field stud- 
ies of the effect of the lower cutoff of fiber strength on 
the failure process of fiber bundle models 0] by con- 
tinuously varying the range of load sharing between the 
limiting cases of completely global load sharing and the 
very localized one [nj • We show that at any range of in- 
teraction there exists a critical value of the cutoff strength 
above which the global response of the bundle becomes 
perfectly brittle. We demonstrate that the crossover of 
the avalanche size distribution to a power law of an expo- 
nent 3/2, when approaching the critical cutoff strength, 
is independent of the range of interaction. Our results 
support the usage of the crossover phenomenon of burst 
sizes in forecasting techniques of imminent failure. 

We consider a parallel bundle of fibers organized on a 
square lattice of size L x L. The fibers are assumed to 
have linearly elastic behavior with identical Young mod- 
ulus £ up to a randomly distributed breaking thresh- 
old. For simplicity, the failure thresholds a t h are as- 
sumed to have a uniform distribution between a lower 
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FIG. 1: a) Failure stress cr c and strain e c of the fiber bundle 
model at zero cutoff el — compared to the critical cutoff 
e l as a function of 7. b) Constitutive curves for GLS 7 = 
(lines) compared to the case of 7 = 5 (triangles) at different 
values of el- 



cutoff strength ol and one with the probability density 
function p(a t h) 
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for a L < <J t h < 1 



0, 



otherwise. 
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Under an increasing external load the fibers break when 
the load on them exceeds the local threshold value <j\ hl 
where i = 1 , . . . , N and N — L 2 denotes the number of 
fibers. Due to the linearly elastic behavior, the failure 
thresholds a\ h can also be expressed in terms of defor- 
mation e\ h = o\ h jE with the cutoff strength £l = &l/E. 
After a failure event, the remaining intact fibers have to 
take over the load of the failed one. In order to give a re- 
alistic description of the load redistribution in FBMs, we 
recently introduced a load transfer function of the form 
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where a a dd denotes the additional load fiber i receives 
after the breaking of fiber j The load increment 

&add decreases as a power law of the distance from 
the failed fiber, where the exponent 7 is considered to 
be a free parameter of the model. The exponent 7 can 
take any values between and 00 controlling the effective 
range of load redistribution between the limiting cases of 
completely global 7 = and completely localized load 
redistribution 7 — > 00 

Under perfectly global load sharing 7 = the macro- 
scopic constitutive equation of the system can be cast in 
a simple analytic form 
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for a L <Ee< 1, 
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where in the following, the value of the Young modulus 
of fibers will be set to unity E = 1. The constitutive 
behavior Eq. (Q of the bundle is perfectly linear up to the 
deformation el since no fibers break in this regime (see 
also Fig. ^) . Due to the breaking of fibers above el, the 
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FIG. 2: Characteristic quantities of the failure process of 
FBM varying the effective range of load sharing 7 and the 
cutoff value of failure strength el- a) critical deformation 
e c , b) critical stress <r c , c) |e c — a c \, d) the average size of 
the largest avalanche (A max ) as a function of el for various 
values of 7. 



constitutive curve a(e) becomes non- linear and develops 
a maximum whose value <J^ LS and position e^ LS define 
the failure stress and strain of the bundle, respectively. 
It follows from Eq. Q that the critical strain is constant 
e^ LS — 1/2 and does not depend on the cutoff strength 
£l , while (Jc increases due to the missing weak fibers 
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Increasing the external load quasi-statically, the break- 
ing fibers trigger avalanches of failure events which cither 
stop after a finite fraction of fibers failed, or become un- 
stable and destroy the entire system. As a consequence, 
the cutoff strength e l can take meaningful values in the 
interval < < e^ LS , since for e^ > ^t LS the breaking 
of the first weakest fiber results in an immediate catas- 
trophic failure of the bundle. 

We explore the effect of the finite cutoff strength el 
on the failure process of FBMs with short ranged load 
sharing by means of computer simulations, redistribut- 
ing the load of broken fibers according to the load trans- 
fer function Eq. Stress controlled simulations were 
carried out on a square lattice of size L = 257 with peri- 
odic boundary conditions varying the cutoff strength el 
of the disorder distribution Eq. JJJ in the interval [0, 0.5] 
at several different values of the effective range of inter- 
action 7 between and 11. To characterize the failure 
process of the bundle at the macro and micro level, we 
determined the critical stress a c and strain e c , the dis- 
tribution D of avalanche sizes A, the average avalanche 
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size (A), and the average value of the largest avalanche 
(A max ). For clarity, we first characterize the behavior of 
the system in the specific case of zero cutoff El = by 
studying the critical stress a c and strain e c of the bundle 
as a function of 7, see Fig. Based on the numeri- 
cal results, three regimes of the failure of FBM can be 
distinguished in Fig. \T£i depending on the range of load 
sharing: for 7 < 2 the range of interaction is infinite in 
the two dimensional embedding space, hence both a c and 
e c take their GLS values er c = 0.25 and e c — 0.5 indepen- 
dent of 7 (see Eq. Q at El, = 0). Increasing the value 
of 7 > 2 the effective range of interaction gradually de- 
creases which lowers the macroscopic strength e c and a c 
of the bundle. In the limiting case of 7 — > 00 the model 
recovers the very localized load sharing, where e c and a c 
take again constant values. According to the numerical 
results, the perfectly localized limit is practically reached 
for 7 > 6, so that in the interval 2 < 7 < 6 a transition 
occurs between the completely global and completely lo- 
cal behavior Fig. ^> demonstrates that for 7 < 2 
(i.e. GLS) the macroscopic failure of the bundle is pre- 
ceded by a strong non-linearity of the constitutive curve 
17(e). At any wider range of load sharing, 7 > 2, the 17(e) 
curves follow the GLS solution Eq. 10, but with lower 
strength values which implies a more brittle macroscopic 
response for short ranged interactions. 

Varying the cutoff strength el at different values of 
7, it can be seen in Figs. Eh, b) that in the long range 
regime 7 < 2 both e c and a c agree well with the analytic 
predictions Eq. Q, i.e. e c = 1/2 is constant while a c in- 
creases with increasing cutoff El- When the load sharing 
becomes short ranged 7 > 2, the increasing macroscopic 
brittleness has the consequence that the curves of e c (el) 
and o c (el) shift downwards as 7 increases and tend to 
a limit curve when the interaction becomes completely 
localized for 7 > 6. It is interesting to note that for 
short range interaction of fibers 7 > 2, not only the fail- 
ure stress (7 C but also the failure strain e c is an increas- 
ing function of El- It is important to emphasize that at 
each 7 there exists a critical value of the cutoff strength 
e c l < sf LS where the failure stress a c and strain e c of 
the system becomes equal to the cutoff strength, i.e. at 
e c l holds e c (e c l ) = u c (e l ) = e c l . At this point the macro- 
scopic response of the bundle becomes perfectly brittle, 
i.e. under gradual loading of the system the macroscopic 
constitutive behavior is linear up to cr c , where the break- 
ing of the weakest fiber gives rise to the collapse of the en- 
tire system (Fig. Eh) . This transition is better illustrated 
by Fig. where the difference 5 = \e c (el) — <t c (el)\ is 
plotted versus El- It can be observed that 8 monotoni- 
cally decreases and becomes practically zero at e c l of the 
given 7. Since the absence of weak fibers gives rise to 
a higher macroscopic strength, the value of £ C L is larger 
than the strength of the bundle e c and a c at zero cutoff 
(see Pig. Eh). 

On the micro level, the failure process is characterized 
by the bursts of fiber breakings, which also show an in- 
teresting behavior when the range of interaction 7 and 
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FIG. 3: Distribution of burst sizes D(A) varying el at dif- 
ferent values of 7: o) 2.0, b) 2.5, c) 3.0, d) 6.0. Crossover 
behavior of D(A) can be observed as El approaches the crit- 
ical cutoff value e l (7) . 



the lower cutoff el are varied. In the GLS regime 7 < 2, 
our computer simulations perfectly recover the analyti- 
cal and numerical results of Refs. 0, E3 (see Fig. Eh): 
for el — the size distribution of bursts D(A) follows a 
power law 

D(A)~A~ a , (5) 

with an exponent a — 5/2. Increasing the value of 
the cutoff el, for small avalanches a crossover occurs 
to a power law of a lower exponent a = 3/2, while for 
large avalanches the original power law with a = 5/2 
is retained. The crossover avalanche size A c increases 
with increasing el and the distribution becomes a sin- 
gle power law of exponent a — 3/2 as El approaches 
ejf is . The crossover to a lower value of the exponent in- 
dicates that due to the missing weak fibers the fraction of 
small avalanches decreases compared to the larger ones. 
This argument is further supported by Fig.|2h and Fig.0] 
which demonstrate that for 7 < 2 both the average size of 
the largest avalanche (A max ) and the average avalanche 
size (A) are monotonically increasing functions of the 
cutoff el- However, when the load sharing gets short 
ranged 7 > 2, both {A max ) and (A) have a maximum at 
the critical cutoff strength. Note that the non-vanishing 
avalanche size above e c l arises due to the strength fluctua- 
tions of the finite bundle so that above e c l the bundle may 
survive a small number of avalanches instead of collapsing 
after the breaking of the weakest fiber. It is interesting 
to note that contrary to GLS, in the transition regime 
2 < 7 < 6, the avalanche size distribution does not show 
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FIG. 4: The average burst size (A) as a function of el for 
different values of 7. In the short range regime 7 > 2, the av- 
erage avalanche size (A) has a clear maximum which coincides 
with the critical cutoff strength e c L . 



a power law behavior for small cutoffs El ~ 0, however, 
when el, approaches the critical value ££(7), the distribu- 
tion of burst sizes D(A) tends again to a power law of an 
exponent a — 3/2 (Fig.©, c). For very localized interac- 
tions 7 > 6 an apparent power law of D(A) is restored 
for el ~ with a relatively high exponent a ~ 9/2, in 
agreement with Ref. (Fig.©). The main outcome of 
computer simulations is that the crossover behavior of 
D(A) to the universal power law -D(A) ~ A -3 / 2 prevails 
at any value of the range of interaction 7 for the limiting 
case of £l — > £l(7)> independently of the original form of 
D(A) at zero cutoff el = (see Fig.©. In spite of the 
relatively large system size L, for short range interaction 
of fibers and El — > ££,(7) the statistics of avalanche sizes 
is rather poor for large avalanches which hinders us to 
make a definite conclusion on the shape of D(A) in this 
A regime. 

An interesting experimental realization of the crossover 



for crackling noise was very recently found in the magni- 
tude distribution of earthquakes in Japan [T^ . Analyzing 
the local magnitude distribution of earthquakes preced- 
ing main shocks, a significant decrease of the Gutenberg- 
Richter exponent was obtained when the lower bound of 
the time window of the analysis is shifted towards the 
catastrophic event 13]. Fracture of ferromagnetic ma- 
terials is accompanied by changes of the magnetic flux, 
which can be recorded as magnetic noise and provides 
information on the dynamics of crack propagation |12| . 
The amplitude, area and energy of magnetic emission sig- 
nals have recently been found to have power law distribu- 
tions with exponents depending on the type of fracture, 
i.e. ductile failure where stable crack propagation occurs 
in a large number of elementary steps is characterized by 
significantly higher exponents than brittle failure, where 
the crack propagates in an unstable catastrophic manner 
breaking the specimen in a few large jumps [12j. Our nu- 
merical results suggest that the reduction of non-linearity 
of the macroscopic response of materials preceding global 
failure when going from ductile and quasi-brittle to brit- 
tle fracture is responsible for the lowering of the crackling 
noise exponents on the micro-level. 

In summary, we carried out computer simulations of 
the failure process of a bundle of fibers with a finite cutoff 
of the fibers' strength, continuously varying the range 
of interaction between the limiting cases of global and 
local load sharing. We showed that increasing the cutoff 
strength el the macroscopic response of the fiber bundle 
becomes perfectly brittle when el approaches a critical 
value £^(7), depending on the range of interaction 7. 
Our numerical results demonstrate the robustness of the 
crossover of the avalanche size distribution D(A) to a 
universal power law of exponent 3/2, irrespective of the 
range of interaction between the material elements. 
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